A noncommutative enumeration problem 



Maria Simonetta Bernabei and Horst Thaler 

Department of Mathematics and Informatics, University of Camerino, 
Via Madonna delle Carceri 9, 1-62032, Camerino (MC), Italy; 

E-mail: simona.bernabei@unicam.it, horst.thaler@unicam.it 



Abstract 

In this article we tackle the combinatorics of coloured hard-dimer objects. This is achieved 
by identifying coloured hard-dimer configurations with a certain class of rooted trees that 
allow for an algebraic treatment in terms of noncommutative formal power series. 
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1 Introduction 

The aim is to count coloured hard-dimer configurations which are objects as shown in 
Figure [1] 



Figure 1: A coloured hard-dimer configuration {n=12) with 2 blue dimers, 1 red dimer 
and 3 inner vertices. 

More precisely, we define a coloured hard-dimer configuration (CHDC) to be a finite 
sequence cr^, of n blue and red vertices, together with coloured dimers on cr^ which 
must not intersect. Here a coloured dimer is an edge connecting two nearest vertices of 
the same colour, see Figure [H The dimer's colour is given by the colour of its boundary 
vertices. In graph-theoretic language CHDCs are a subclass of labeled graphs whose 
vertices and edges carry one of two possible labels. Without loss of generality we may 
assume that a CHDC is a subset of R and its vertices belong to Z. We include also the 
empty CHDC, i.e. the configuration when no dimers are present. 

The enumeration problem we have in mind goes as follows: Given a sequence an 
and a triple of nonnegative integers {i,j,k), where i,j count the numbers of blue, red 
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dimers and k the number of inner vertices, i.e. vertices that are connected to their left 
and right vertex. How many CHDCs of a given type {i,j, k) are there on an, where an 
as well n may he arbitrary ? 

The motivation for studying this kind of objects is due to their appearance in 
causally triangulated (2+l)-dimensional gravity. It was shown in [T], by using spe- 
cial triangulations of spacetime, that the discrete Laplace transform of the one-step 
propagator can be expressed as follows 

Z{u, f , w) = ^ e~'^"Z„(M, f , w), (1) 

nGN 

with 

Zniu, v,w) = Y u ^ T and Zl'^Hu, vM= V 

Cn ^ ' CHDCs D|<,„ 

In reference [T] it is explained how the parameters ^^u^v^w are related to physical and 
geometrical constants. The exponents |-D|fe, |-D|r; count the number of blue dimers, 
red dimers and inner vertices, respectively. Now, if we want to find explicit expressions 
for the function Z above, we need to know how many CHDCs, for a given sequence cxn, 
give rise to the same triple of integers (|-D|b, |-D|r, |nD|). This is exactly what we are 
after in this article and we shall find an algebraic solution to the enumeration problem 
as stated above. The strategy is to identify the set of CHDCs with a certain class 
of trees the elements of which can be encoded into monomials of a noncommutative 
formal power series ^^'^^x, where the x's are words over the alphabet {6, r}. This 
means that the indeterminates h and r are supposed to be noncommutative and a 
particular word x just corresponds to a particular sequence of blue and red vertices an- 
The coefficients Ox collect the information of CHDCs of a given type on o"„. 

In section [2] we show how to make CHDCs into monomials and thereby solve the 
enumeration problem from above. In section |3] we give an alternative way to express the 
coefficients Ox by means of a linear representation of the indeterminates 6, r. The formal 
power series that allow such a representation are called recognizable. This notion is 
explained in Appendix A and in Appendix B the matrices defining the representation 
are given. Finally, we use this result to investigate the growth behavior of the number 
of CHDCs. 

2 From CHDCs to noncommutative series 

In this section we shall establish a bijection between CHDCs and a certain class of 
trees. This is achieved by first identifying CHDCs with a particular class of graphs 
denoted Gshcd which in turn will be identified with the corresponding class of trees 
denoted Th^d- 

Bisection CHDCs <^=^ G/jc^ : For a given CHDC we endow its vertices with the order 
inherited from Z. An extra vertex x is added in the first place, which becomes the root 
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(a) 




(b) 



Figure 2: Illustration of the one-to-one correspondence between a coloured hard-dimer, 
its graph and its tree. 



of the graph. Then all the vertices which are next neighbors are connected through an 
edge. Finally, the vertices which in the original CHDC are the starting and end points 
of a dimer are linked via an extra edge to the right, as seen from the root (see Figure 
El^a)). We note that only the closest vertices of the same colour may get an extra link 
and that these become trivalent, unless we consider the last vertex, in which case it 
becomes bivalent. The set of rooted graphs we obtain in this manner is denoted Ghcd 
and it is obvious that we have established a bijection between the latter and the set of 



Bijection Ghcd -v^ T^hcd '■ When a graph in Ghcd is given, we move beginning from 
the root in clockwise direction and cut an edge if the graph remains connected. The 
cut edge is replaced by two edges ending in univalent vertices called buds (black arrow) 
and leafs (white arrow), respectively. The procedure is repeated until we reach the last 
vertex (see Figure [2]^b)). We note that buds and three vertices are always connected by 
one edge only whereas leafs and three vertices may be interlaced by a certain number of 
bivalent red or blue vertices. The set of rooted trees we obtain following this procedure 
is denoted Thcd- Viceversa, when a tree in Thcd is given we may, starting from the root, 
move in clockwise direction and merge the bud leave pairs to edges getting a graph in 



We have thus established a bijection between the set of CHDCs and the set Thcd- 
The next goal is to find an algebraic formulation of the enumeration problem for ele- 
ments in Thcd- For this we adapt ideas from [3] and assign the charges +1,-1 and 
to leaves, buds and univalent coloured vertices, respectively. With this choice any tree 
has charge and any subtree different from a bud has charge or -t-1. For this reason 



CHDCs. 
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we distinguish between S (respectively, R) trees, iff: 

(a) their total charge is (respectively, 1); 

(b) any descendant subtree different from a bud has charge or 1. 

Furthermore, we subdivide these two groups into Sb, and S*,., Rr trees depending on 
whether the first vertex after the root is blue or red. 

Now we are able to write down a system of equations for the four types of trees which 
will turn out to be linear. Since we are interested in counting the number of CHDCs 
on a fixed sequence of blue and red vertices we assign noncommutative variables h 
and r to blue and red vertices, respectively. Furthermore, for a given sequence an we 
group together all CHDCs with the same numbers |-D|b, |-D|r and | nD|. This can be 
achieved by encoding a blue or red dimer through the assignment of commutative vari- 
ables 63 and r3, respectively. Particularly, the variables 63 and are assigned to those 
trivalent blue and red vertices, respectively, which possess a bud leg pair as subtree. 
Finally, to every bivalent vertex sitting in between a three vertex and a leaf we assign 
a commutative variable y. 

We are thus naturally led to consider the trees Sb, Sr,Rb, Rr as elements of the 
noncommutative algebra K{{b,r)) of formal power series where K is the commutative 
polynomial ring K = F[63,r3,?/] (with F being any subfield of M which contains Z). 
For example, the tree and the corresponding graph in Figure [2] read 

b'^r^y^ r*b*r*r*b*r*b*b*r*b*r*b, 

where * stands for the noncommutative product in K{{b,r)). As we can read off from 
Figure 3 our trees have to satisfy the following consistency equations: 

Sb = b + b*{Sb + Sr) + bsb*Rb 

Sr = r + r * {Sb + Sr) + r^r * Rr 

R^ = })-^b*{Sb + Sr)+yr*{l-yr)~'^*b*{Sb + Sr) + yr*{l-yr)~'^*b 

R^ = r + r * {Sb + Sr) +yb* {1 -yby^ *r * {Sb + Sr) +yb* {1 -yby^ *r 

Substituting the expressions of Rb and Rr into the equations for Sb and Sr we get by 
elementary algebra the following system of equations 



Sb= [6 + 63 6^ + 63?/ 6 * r * (1 - yr)'^ * b] *Sb + Ab * Sr + 



b 



Sr = [r + r^r"^ + r^y r * b * {1 - yb)'^ * r] *Sr + Ar * Sb + Ar, 



that is 



(1 - Ab) *Sb-Ab* Sr = Ab 



(2) 
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The system ([2]) can be solved easily with the result 



= (1-A^*(1-A,-Ah)-^-1. 



(3) 



By symmetry we also get 



Sr = (1 - A) * (1 - A, - - 1. 
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Figure 3: The possible building blocks from the perspective of an arbitrary vertex. The 
corresponding building blocks for Rr are obtained by exchanging colors and variables 
(6 <(-)■ r). 
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Proposition 2.1 Let Sf, = Xlx'^x^ be the formal series given as part of the unique 
solution of Ij^ in K{{b,r)). The sum runs over all words x G {b,r}* starting with 
letter b and the coefficients = {St, x) are given as finite sums 

Ox = ^ m^j- fc(x) b^r^^y'', 

jJ,fcGN„ 

where the multiplicities rriij^ki'^) count the number of coloured hard-dimer configura- 
tions D on ^ for which (|Z)|b, | H D|) = {i,j,k). In particular, the evaluation 
ax(&3 = l,r3 = 1, w = 1) gives the number of CHDCs on the sequence (T„ corresponding 
to X. The symmetric assertion holds for Sr- 

Proof. Let # be one of the subscripts b, r and decompose S# = where 5*^ 

comprises the sum of terms in with |x| = In the same manner we write = 
-^#- elementary computation shows that 5*^ = {higher order terms} and 
hkewise for R^, which implies that a solution to ([2]) gives rise to sequences (S'^)„>_i 
and {R^)n>-i obeying the following recursive equations {n > 1): 

= b + b*is;^-^ + si!-^) + bsb*R]:-\ 

S;^ = r + r*{SJ^-' + S';-') + r;r*R'^-\ 

n~2 

Ri = b+b* {sj^-^ + ^r') + XI * ^ * (-^r'"^ + -^r'"^) 

fc=i 

n-2 

R^ = r + r* {Sl'-^ + S;'-') + X yH'^ * r * + S;'-'-^) 

k=l 

+ *r, (4) 

with initial conditions 5*^ = = for c G { — 1,0} and = R^ = Vice 

versa, for sequences (S'^)„>_i, (-RJ)„>_i obeying (jl]) and related initial conditions, 
their sums S'# := Yl'^=-i ^# ^'^^ ^# ■= J2'^=-i will satisfy (j2]). Therefore we may 
equivalently look at the system (jl]). But the elements contain precisely the sum of 
those terms from K{{b,r)) which are algebraic counterparts of trees in Thcd that have 
n coloured vertices the first one of which is 7^. This is due to the fact that there is 
a one-to-one correspondence between trees from Thcd and trees of charge that are 
constructed recursively according to the building blocks from Figure 3. For a fixed x, 
the coefficient Ox is a sum of elements in K which stem from trees that are constructable 
for the particular sequence of blue and red vertices given by x. Since the indeterminates 
63, r3,?/ are commutative, trees will contribute the same term b^r^y^, whenever the 
numbers of vertices corresponding to a b-^,r-^,y indeterminate, coincide. Moreover, a 
glance at the bijection CHDCs T^cd shows that {i,j,k) = (|-D|b, \D\r, |ni^|), which 
proves the assertion. □ 

As immediate consequence we get 
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Corollary 2.1 Let Sf, = X^x'^^x be the formal power series representing the solution 
of The discrete Laplace transform from (QP can be expressed, using Proposition 
\2.1\ by means of 

Z^f{u, V, w) = [(ax(&3 = u,r3 = v,y = w) + a^ib^ = v,r^ = u,y = w)]. 

In the latter equality we have evaluated the coefficients Ox at points u,v,w G C and 
have exploited the symmetric relation between Sf, and Sr- 

3 The solution as a recognizable series 

In order to find the coefficients in we might start directly from ([3]). This, however, 
necessitates to invert and multiply noncommutative power series. A recursive proce- 
dure or the use of Cauchy's product would mean an exponentially growing number 
of operations. Therefore, if we'd like to employ numerical or symbolical programs on 
computers, it is worthwhile making some more effort to find an explicit representation 
for the coefficients Ox of Sb by means of matrices. This is in fact possible as Sb turns 
out to be a recognizable series, see Appendix A. 

Proposition 3.1 The series Sb is a recognizable series with representation 



given by = Bi, /ii(r) = Ri, where the matrices Bi, Ri and tuples Ai, 71 G K^^, are 
given in Appendix B. An analogous statement holds for Sr with ^2(1^) = -Bi(&3 ^ fs) 
and ^2{b) = Ri{rz ^ b^). Since the property of being recognizable is preserved under 
algebraic operations we also have that Sb + 5*^ is recognizable. 

Proof. By Proposition lA. II we have to show that Sb is contained in a finitely generated 
(f.g.) stable i^-submodule of K{{b,r)). Hence we shall verify one by one that the 
following elements are recognizable series: Ar, 1 — Ar — Ab and Sb- Two basic rules are 
needed for proving the statement. For a letter a and any two formal series P, Q we 
have: 



fii:{b,ry^K 



19x19 



a-\PQ) = {a-^P)Q + (P, l)(a-^Q) 



and for a proper series Q 



a-\Q*) = {a-^Q)Q*. 



So ffist let's have a look at 




It is not difficult to see that the finitely generated (f.g.) iC-submodule 



= Kl® KS2 ®Kr® KS2 *r ® Kb® Kb* S2® Kb* S2*r, 
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where 5*2 = (1 — yb) ^, is stable and contains the term b * {1 — yb) ^ * r. In the same 
vein one shows that the i^-submodule U2 generated by 

Si = 1, 5*2 = (1 - yby^, S3 = r, S4 = S2 * r, S5 = b, Sq = b * S2, S7 = b* S2*r, 

Sg = r * b, Sg = r * b * S2, Siq = r * b * S2 * r, Su = r^. 

is stable, f.g. and contains = S3 + rsSii + rs?/ Siq. 
Enlarging U2 by the following basis elements 

5*12 = (1 - yry^, Si3 = S12 *b, Si4 = r* Su, S15 = r * S12 * 6, 5*16 = 6 * r, 

Si7 = b* r * S12, Sis = b* r * S12 * b, Siq = 6^ 
we get a f.g. stable i^'-submodule U3 which contains 1 — Ar and 

1 - Ar - Ab = 1 - S3 - r3 Sii - r3y Sw - 5*5 - 63 Sig - b3y Sis- 

Finally, we define P := Ar + Ab, P* := {Ar + Ab)* and the i^'-submodule generated 
by the basis elements 

Ti := 1, T2 = P*, T3 := 52 * P*, T4 := 53 * P*, . . . , T20 := Si, * P* . 

An elementary, yet tedious calculation shows again that f/4 is a f.g. stable submodule 
containing 

Sb={l- Ar) * (1 - Ar - Ab)-' - 1 = (1 - Ar) * P* - 1 

= P*-S3*P*- r3Sii * P* - r3ySio * P* - 1. 
The corresponding statement for Sr follows once more from symmetry considerations. 
□ 

The result above allows to study the asymptotic behavior of the number of CHDCs 
on cT„ when n tends to infinity. According to Proposition (13.1 p Sb + 5*^ is a recognizable 
series and the evaluation at 63 = rs = ?/ = 1 again gives a recognizable series. This 
entails that for 

M = Sbib3 = l,r3 = 1,1/ = 1) + 5,(63 = l,r3 = 1,2/ = 1) 

there exists a representation fi' : {b,r}* — t- M.'^^'^ and vectors A', 7' G such that 
(M, x) = A'"^ • /i'(x) ■ 7'. This representation is given by 

/x'(6) = B{b3 = l,y = l) and /i'(r) = R{r3 = l,y = l) 

and the vectors by 

A' = A(63 = l,r3 = l,2/=l), 7' = 7, 
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see Appendix B. Note that (M, x) counts the number of CHDCs on the sequence of 
blue and red sites corresponding to x, irrespective of whether this sequence starts with 
b or r. 

It is convenient to discuss the asymptotics of CHDCs within the framework of 
ergodic dynamical systems. To be precise, let B' := fi'{b),R' := /i'(r) and define u 
as the probability measure on {B',R'} given by z^(-B') = i^iR') = |- Then as ergodic 
dynamical system J^, P, 9) we choose 

• Q = {B', R'}^, i.e. the set of sequences u = {ui, U2, ■ ■ .) with Ui G {B', R'}. 

If gn{oj) = Un denotes the coordinate maps then should be the cr- algebra gen- 
erated by the gnS. 

• P is the product measure P = H^i ^• 

• 6^ is the shift operator 9{ijj) = {ijJ2, us, . . .). 

Upon turning {6, r}^ into a probability space that is isomorphic to (fi, J-", P) we may 
state the following theorem 

Theorem 3.1 For u G {b, r}^ let Un '■= (cui, . . . , oon) o,nd let Dn{uj) denote the number 
of CHDCs on ujn- There exists a finite constant random variable a on ({6, r}^, J^, P) 
such that for almost all u G {b, r}^ 

— In Dn{uj) — )■ a, as n —)■ 00. 

n 

This means that the number of CHDCs grows exponentially and that asymptotically 
this growth rate is the same for almost all uj. 

Proof. If we define the matrix-valued random variables by 

Xn{uj) = gi{u)g2iuj) ■ ■ ■fi'„(a;), 

then 

F„ = -ln(A'^-X„-y), (5) 
defines a subadditive sequence of random variables, i.e. 

^m+n < y-rn ° 0'^ + Yn, for all m and n. 

To see subadditivity we first note that with the natural identifications a;„ ^ x : 
(wi, . . . , Un) = wi * ■ ■ ■ * = X we may write A'"*" ■ ■ 7' = (M, x). Writing x = 

Xm * Xm', with |xm| = TTL and |xm'| = TTi' , subadditivity follows from the combinatorial 
fact that (M, x) > (M, Xm)(M, x^')- Now Kingman's subadditive ergodic theorem, see 
[6] and [H Th. IV.1.2], applied to our ergodic dynamical system {fl,J^,F,9) implies 
that there exists a random variable /3 such that 

lim ( — ) = f3 = E[f3], P-a.s.. 

n-j-oo \ n J 



10 



In addition one has /3 = inf„ (E[n ""^Vn])- A priori it is not clear whether /3 > — oo 
which we clarify now. By Jensen's inequality we find 

E [n-'^Yn] > -n-^ ln(A'^ ■ E[X„] ■ 7') = -n'^ ln(A'^ ■ ■ 7'), (6) 

B' -\- R' 

where Z = . From the very definitions we deduce easily that 

A'^ ■ . 7' = 2A'^ ■ ^B,{h = l,y = l) + RM = l,y = l)'^ \ ^, 

An explicit calculation shows that the dominant eigenvalue of the nonnegative matrix 

^ ^1(63 = 1, y = 1) + Ri{n = 1, y = 1) 

2 

is 3/2 with algebraic multiplicity one. Using an appropriate representation for the ma- 
trix powers like that employed in [51 Th. 8, Ch. Ill] shows that lim„_^oo —n~^ In /(ra) = 
3 

— -. This in turn implies E [n^^F^] > const., Vra, for some n-independent finite con- 
stant and consequently /3 = inf„(E[n^-'^l^]) > const. We have thus established almost 
sure convergence of 

— In Dn = ^ —/3 =: a, as n — )> 00. □ 

n n 

A Formal power series 

In this Appendix we recall some basics related to recognizable power series. For more 
details the reader is referred to [2] or [7]. Suppose we are given a semiring K and an 
alphabet A with its associated free monoid A* and unit 1. A formal noncommutative 
power series S or simply a formal series is a map 

S:A*^K. 

The image by S* of a word x is denoted by (5*, x) and is called coefficient of x in S. The 
coefficient (5, 1) is called constant term and a formal series S whose constant term is 
zero is called proper. For x = xi * ■ ■ ■ * a;„ we define |x| = n to denote the length of 
the word x with the convention |1| =0. The set of formal series with coefficients in K 
is denoted by K{{A)). We shall consider K{{A)) as being equipped with the natural 
semiring structure where addition and product are given by 

(5 + T)(x):=(5,x) + (T,x), (5T)(x):= (S,y){T,z). 

y, zGj4* : y*z=x 

Moreover, there are left and right external operations of K on K{{A)) which are given 
by 

(/cS*, x) := fc(S', x); {Sk,x.) := {S,x.)k. 
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The subalgebra K{A) is defined to comprise those formal series with only finitely many 
terms different from zero. It is obvious that with the operations defined above K{{A)) 
becomes an algebra when if is a ring. The semiring K{{A)) can easily be turned into 
a complete topological semiring. This is achieved by introducing the following metric. 
First, let 

K : K{{A)) X K{{A)) ^ No U cx), 
k{S, T) := inf {n G No | 3 x G A* : |x| = n and (S, x) ^ (T, x)} 
and then define the metric 

d:K{{A)) X K{{A)) ^R+, 
d{S,T) := • 

It is easy to see that with respect to this topology a sequence {Sn)neN converges to 5" if 
and only if Vm G No 3 no G N : [Vn > no and Vx with |x| < m : (S'„, x) = {S, x)] . This 
fact then implies that every Cauchy sequence is convergent and it becomes natural to 
write elements S G K{{A)) in the form 

S = (^,x)x. 

For a proper series S it is possible to define the star operation hj S* := X]^o 

Definition A.l (a) The rational operations in K{{A)) are the sum, product, the two 
external operations and the star operation. A subset of K{{A)) is rationally closed if 
it is closed with respect to the rational operations. The rational closure of a subset E 
of K{{A)) is the intersection of all rationally closed subsets of K{{A)) which contain 
E. When K is a ring the star operation and inversion play the same role since S* = 

ii-sr\ 

(b) A formal series is called rational if it is contained in the rational closure of K{A). 

(c) A formal series S is called recognizable if there is some n > 1, a homomorphism of 
semirings (or simply a representation) 

/i : A* ^ if"^", 

where if"-^" carries its multiplicative structure and tuples A,7 G if", such that for all 
words X 

(5,x) = AT-Mx)-7- 
A i^'-linear operation of A* on K{{A)) can be defined by setting 

y-i^:= 5^(^,y*x)x, y e A*. 

A subset M of K{{A)) is called stable if for all x G A* we have x-^(M) C M. The 
following proposition is the core for proving Proposition 13.11 
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Proposition A.l A formal series S is recognizable if and only if there exists a stable 
finitely generated left K-submodule of K{{A)) which contains S. 

In addition the following equivalence holds 

Theorem A.l (Schiitzenberger) A formal series is recognizable if and only if it is 
rational. 



B The matrices B and R 

Let (Ti, . . . ,T2o) be the basis elements spanning the submodule as defined in the 
proof of Proposition 13.11 Furthermore, let B =: fii{b),R =: fii{r) be the matrices 
given by b^^Ti = Yl'fLi^ij'^j f~^Ti = Yl'j=i -^ij'^j ■ turns out that the first row 
and first column of B and R contain only zeros which implies that we may reduce the 
representation /ii to the subspace Vi := spanj^{T2, . . . ,T2o} where it is determined by 
fii{b) = Bi,fii{r) = Ri as 
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Finally, define A, ^ hy Sb — Ylfli ^i^i and 7^ = (Tj, 1) and set 

Ai := (A2, A3, . . . , A20), 7i := (72,73, • • • ,720), 

then 

Ai - (10 -1 000000 -rsy -rg ) , 
71 = (1100000000010000000). 



The corresponding representation for Sr can be found easily by exchanging the 
letters b and r. For this we define the new basis elements 

fi := T,{b ^ r)2<i<2o, 

which are linearly independent of the Tj's, and set V2 := span;^{T2, . . . , T20}. 
Furthermore we put A2 := Ai(r3 i-)- 63) and 72 := 71- Thus the series Sr is recognizable 
with representation 

//2:{6,rr^X^9x^ 

given by 

Ii2{b) = B2 := i?i(r3 ^ 63), /i2(r) = R2 := Bi(&3 ^ ^3), 

and tuples A2,72 G K^'^ . It is readily seen that also + Sr is recognizable with its 
representation defined on K ® V2: 

//:{6,rr ^i^38x38 

is specified by 



^{h) = B = [ ] and /i(r) = i? 



Ri 

i?2 



and the tuples A, 7 e K^^ © K^^ are given as A = (Ai, A2) and 7 = (71, 72). 
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